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ABSTRACT 



Aims. We investigate the mineralogy and dust processing in the circumbinary discs of binary post-AGB stars using high-resolution TIMMI2 
and SPITZER infrared spectra. 

Methods. We perform a full spectral fitting to the infrared spectra using the most recent opacities of amorphous and crystalline dust species. 
This allows for the identification of the carriers of the different emission bands. Our fits also constrain the physical properties of different dust 
species and grain sizes responsible for the observed emission features. 

Results. In all stars the dust is oxygen-rich: amorphous and crystalline silicate dust species prevail and no features of a carbon-rich component 
can be found, the exception being EPLyr, where a mixed chemistry of both oxygen- and carbon-rich species is found. Our full spectral fitting 
indicates a high degree of dust grain processing. The mineralogy of our sample stars shows that the dust is constituted of irregularly shaped 
and relatively large grains, with typical grain sizes larger than 2 /jm. The spectra of nearly all stars show a high degree of crystallinity, where 
magnesium-rich end members of olivine and pyroxene silicates dominate. Other dust features of e.g. silica or alumina are not present at 
detectable levels. Temperature estimates from our fitting routine show that a significant fraction of grains must be cool, significantly cooler 
than the glass temperature. This shows that radial mixing is very efficient is these discs and/or indicates different thermal conditions at grain 
formation. Our results show that strong grain processing is not limited to young stellar objects and that the physical processes occurring in the 
discs are very similar to those in protoplanetary discs. 
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1. Introduction 

Post-AGB stars are low- and intermediate-mass stars that 
evolve rapidly from the Asymptotic Giant Branch (AGB) to- 
wards the Planetery Nebula (PNe) phase, before cooling down 
as white dwarves. During the previous AGB phase, the star 
has undergone severe mass loss, leaving behind a slowly ex- 
panding circumstellar dust shell. Depending on the dilution 
in the line-of-sight, the expected spectral energy distribution 
(SED) is double-peaked, with the first UV-visible peak indica- 
tive of the central star and the second infrared peak coming 



* Based on observations obtained at the European Southern 
Observatory (ESO), La Silla, observing program 072.D-0263, on ob- 
servations made with the 1.2 m Flemish Mercator telescope at Roque 
de los Muchachos, Spain, the 1.2 m Swiss Euler telescope at La 
Silla, Chile and on observations made with the SPITZER Space 
Telescope (program id 3274), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract with 
NASA. 



from thermal emission of the cool dust in the circumstellar 
environment (CE). The release of the IRAS All-Sky Survey 
mission made a large-scale identification o f post-AGB stars on 
the basis of their infrared colours possible jLlovd Evaiislll985 



Hrivnak et al.lll989t lOudmaiier et al. 1992^. More recent sur- 
veys on post-AGB stars include Tszczerba et al. (2007), who 
present the latest catalogue of Galactic post-AGB stars with 
almost 400 objects. 

IRAS colour-colour diagrams revealed a large sample of 
stars that did not show the expected d ouble-peaked SED 
(lLloydEvanslll985 : iDe Ruvter et al ] l2006h . ' nstead, they show 
a strong near-IR excess, pointing to the presence of hot dust 
in the system, while the central stars are cuiTently too hot 
to have an ongoing dusty mass loss. Correlations with opti- 
cal databases of variable stars point to the presence of pul- 
sating stars that also show this specific IR-excess. These 
RVTauri stars are luminous evolved stars that cross the 
Population 11 Cepheid instability strip and pop ulate a well- 
defined part of the IRAS colour-colour diagram (ILlovd Evans 
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19851; iRaveendranll 19891: iLloyd Evansl ll999h. Recent interfer- 
ometric studies ( Deroo et al. 2007allb ) prove that the circum- 
stellar emission originates from a very compact region and 
that the SEDs of these stars are well modelled with a pas- 
sive 2D disc model ( DuUemond & Dominikli2004l: lOeroo et al.l 



2007at iGielen et al.ll2007h . The spectral slope of the submil- 
limetre SED points to the presence of large, up to and cm 
sized , grains in these discs (iDe Ruyter et al iGielen et al.l 

2007 1). which have settled to form a cool disc midplane. The in- 



ner radius of these discs is determined by the dust sublimation 
radius. The hot inner rim is puffed up by gas presssure from the 
central star and radiates mainly in the near-lR, while the outer 
parts of the disc can be strongly flared and are responsible for 
the strong absorption and re-radiation of the stellar light. In a 
few cases the direct confirmation of the Keplerian kinematics 
of the circumstellar discs comes from resolved interferometric 



CO measurements dBujarrabal et al.lll988ll2005ll2007h . Recent 
studies have shown that these post-AGB stars are v ery likely 



all binaries (iM aas e t al]|2002l 12003 



Van Wincke II2OO7I: IGielen et al.ll2007h . 



Van Winckel et al. 1999; 



Famous examples of these objects include HD 44179, 
the central star of the Red Rectangle. It is a carbon- 
ric h post-AGB objec t for w hich the binarity was proven 
by 



post-AGB objec t for w hich the 
Van Winckel e t al. (1995). The star is surrounded by an 



oxygen-rich disc ([Waters et al. 1998) and is resolved in ground 



based high spatial resolu tion imaging (IRoddier et al.l 11995 



iMen'shchikov et al.l 120021) as well as in HST optical images 



(Cohen et al . 20041). The disc is also resolved in interferomet- 
ric CO measu rements, where the Kep lerian velocity of the disc 
was detected dBujarrabal et al ] l2005l) . 

Infrared spectroscopy is the ideal tool to study the 
physico-chemical characteristics of the circumstellar mate- 
rial, since it samples resonances of the most dominant 
dust species and important ro-vibrational bands of dominant 
molecules. Spectroscopic data obtained with the Infrared Space 
Observatory (ISO) allowed for the first time to study the miner- 
alogy of circumstell ar environments o f objects i n differen t evo- 
lutionary stages (e.g. lWaters et al 1IT99 6; Waelkensetaill996). 
More advanced studies of the CE of yo ung stellar objects (e.g. 



van Boe kel et al.ll2005 | ; |Li s s e et al.ll2007i) and mass-losing ob 



jects (Molst eTet alj 2002al b 3) further identified the dominant 
dust species and grain sizes of the circumstellar dust. 



2003 



Infrared spectral studies of post-AGB objects dVan Winckel 
and references therein) allowed for the detection of 



different CE chemistries. The more typical C-rich post-AGB 
stars are characterised by a strong amorphous SiC feature at 
1 1.3 jum. In these star s, a strong unident ified 21 fim feature can 
sometimes be found ( Kwoket alj 19891) . In addition to the de- 
tection of O-rich post- AGB stars, dominated by silicate species, 
surprisingly also mixed chemistries were found, showing fea- 
tures of both O-rich and C-rich dust species. 

Young stellar objects are the natural environment to study 
circumstellar disc physics, but previous studies of the miner- 
alogy of discs a round the few brigh t est infr ared evolved ob- 
jects with ISO (' Molster et aP 1 19991 l2002 a'b.'c) show strong 
dust processing, with oxygen-rich and highly crystalline dust. 
Molster et al. (I2002allbll3) found that a high degree of crystalli- 



sation is indicative for the presence of a stable disc, and not a 
dusty outflow. 

In our pilot study dCielen et al.ll2007l) . based on SPITZER- 
IRS and ISO data, we investigate the mineralogy and spec- 
tral energy distribution of two post-AGB stars, RUCen and 
AC Her The spectra of these stars are extremely similar and 
show a strong resemblance to the infrared spectrum of the 
solar-system comet Hale-Bopp and young stellar objects, such 
as HD 100546. The observed crystalline emission features 
are well modelled using magnesium-rich crystalline silicates. 
The grains are irregularly shaped, with grain sizes of 0. 1 fim 
and 1.5 /vm. Both hot and cool grains are necessary to re- 
produce the observed spectrum. The spectral energy distribu- 
tions of both stars were fitted using a 2D passive disc model 
dDuUemond & Dominik 2004). The discs surrounding these 
objects start at 35 AU from the central star, which is well 
beyond the dust sublimation radius. The outer radius is less 
well constrained, but extends till about 300 AU. The discs are 
strongly flared. 

In this work we study the mineralogy of the discs in a 
large sample of post-AGB binaries and perform a detailed fit- 
ting of the observed emission features. We have observed 21 
binary post-AGB objects with the SPITZER-IRS spectrograph 
and look for relations between the dust parameters and stellar 
characteristics. 

The outline of this paper is as follows: in Sects. |2]and [3] 
we introduce our programme stars and the different observa- 
tion and reduction strategies used. Section |4] contains the con- 
struction of the spectral energy distributions and the total ex- 
tinction determination. In Sect. |5] we give a general overview 
of the spectra and the observed emission features. The profile 
and full spectral fitting is presented in Sects. |6]and [T] The dis- 
cussion and conclusion of our analysis are given in Sects. [8]and 

m 

2. Programme stars 

From the total sample of 5 1 publis hed binary post-AGB stars 
likely surrounded by a stable disc dDe Ruyter et al.ll2006l) . we 
selected the 21 stars with fluxes below the saturation limit of 
the SPITZER-IRS spectrograph. Spectral types range from A 
to M. 

Our radial velocity monitoring program, using the 
CORALIE spectrograph attached to the 1.2 m Swiss Euler 
Telescope at the ESO La Silla Observatory, is still ongoing but 
we have already found orbital parameters for 15 of our 21 pro- 
gram stars. Orbital periods ran ge from 100 up to more than 
2000 days dVan Winckel 120070 . The objects for which we do 
not yet have orbital parameters, have pulsational amplitudes 
which are too large for a straightforward detection of the binary 
motion. The results of the radial velocity monitoring program 
will be subject of a forthcoming paper 

The orbital periods indicate that strong binary interaction 
must have occurred during the late stages on the AGB. The bi- 
naries are now not in contact but the orbits are too small to ac- 
commodate a full-grown AGB star The discs are all circumbi- 
nary since all orbits lie well within the dust sublimation radius. 
The distribution of the mass-functions gives a range in minimal 
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mass for the companion between 1 Mq and 2 Mg. The compan 
ion stars are likely to be unevolved m ain sequence stars (e.g 
Gielen et al l2007HVan Winckelll2007h . 

These binary post-AGB stars are c haracterised by a de 



The spectra were extracted from the SSC raw data 
pipeline version S13.2.0 products, usin g the c2d Interactive 



pletion pattern in their photosph er es (IGiridhar et al 



1998", '2000; 'Gonz alez et aLlll997all5 IVan Winckel etal 



1994 



1998 



Maas et al. 2005). This abundance pattern is the result of gas- 
dust separation followed by r eaccretion of the g as, which 
is poor in refractory elements. I Waters etal (1992) proposed 
that the most likely circumstance for this process to oc- 
cur is when the dust is trapped in a circumstellar disc. 
Photospheric chemical studies of post-AGB candidates in the 
LMC have revealed that there also, de pletion patterns are com- 
mon ( Revniers & Van Winckell 2007b . T he observations in the 
recent release of the SAGE database (" Meixner et al. 2006 ) 
showed that these depleted LMC sources have infrared ex- 
cesses similar to the Galactic binary post-AGB stars, and are 
therefore thought to be disc sources as well. 

The disc formation itself is badly understood. Possible for- 
mati on scenarios of the discs includ e a wind capture scenario 
("e.g. iMastrodemos & MoiTis 1999), or a formation scenario 



through non-conservative mass transfer in an interacting bi- 
nary. In the first scenario, the AGB wind is captured by the 
companion. In the second scenario, which is still not very well 
explored theoretically, the disc formation precedes the dust- 
grain formation and it is likely that the thermal history of the 
grains was very dififerent from that in normal AGB winds. This 
could lead to very different chemico-physical properties during 
formation. The sizes of the orbits suggest this scenario to be 
more likely. 

We may witness the formation of a circumbinary disc by 
Roche-Lobe overflow in the evolved binary system SSLep. 
The optical/IR interferometric observables can be best under- 
stood, assuming a Roche-lobe filling M s tar in a system with a 
circumsystem disc dVerhoelst et al ] |2007h . 

Our program stars are therefore all likely binaries which are 
surrounded by a dusty stable disc. This disc seems to have an 
important impact on the objects and in this paper we focus on 
the analysis of the IR spectra as probes for the disc physics. 



3. Observations and data reduction 

3.1. SPITZER 

High- and low -resolution spectra for 21 post-AG B stars were 



obtained using the Infrared Spectrograph (IRS: iHouck et al 
2004) aboard the SPITZER Space Telescope dWerner et al 
20041) in February 2005. The spectra were observed using com- 



binations of the short-low (SL), short-high (SH) and long-high 
(LH) modules. SL (/l=5.3-14.5jum) spectra have a resolving 
power of R=^/ A ^ ~ 100 , SH (.l=10.0-19.5yum) and LH 
(/l=19.3-37.0//m) spectra have a resolving power of ~ 600. 
Exposure times were chosen to achieve a S/N ratio of around 
400 for the high-resolution modes, which we complemented 
with short exposures in low-resolution mode with a S/N ratio 
around 100, using the first generation of the exposure time cal- 
culator of the call for proposals. 



Analy sis reduction so ftware package dKessler-Silacci et al 



^2006; iLahuis. F. et al] 12006). This data processing includes 
bad-pixel correction, extraction, defringing and order match- 
ing. Individual orders are corrected for offsets, if necessary, by 
small scaling corrections to match the bluer order. 

3.2. TIMMI2 

For stars where we lack the SPITZER IRS-SH observations 
we obtained additional ground-based N-band infrared spectra 
with the Therrnal Infrared Multi Mode Instrument 2 (TIMMI2, 
Reimann et al. I l2000t iKaufl et al.ll2003h . mounted on the 3.6 m 
telescope at the ESO La Silla Observatory. The low-resolution 
{R ~ 160) N band grism was used in combination with a 
1.2 arc sec slit, the pixel scale in the spectroscopic mode of 
TIMMI2 is 0.45 arcsec. For th e reduction of the spectr a we 



used the method described in Ivan Boekel et al.l (120051) . We 



scaled the TIMMI2 spectra to the SPITZER spectra and found 
a very good agreement in spectral shape between the two data 
sets. 



4. Spectral energy distribution 

For all sample stars spectral energy di stributions were update d 
from the photometric data as given by lDe Ruyter et al. I (l2006t) . 
The resulting SEDs are presented in Figs. IB. II and IB. 21 The 
total extinction E(B - V),ot was determined by dereddening 
the observed pho t ometr y, using the average extinction law of 
Savage & MafliisI d 1979b . The relation between E{B - V),ot and 
Av is given by Ay - Ry x E(B - V) , with a typical value 
for Ry = 3.1 dSavage & Mathislll979l) . Minimising the dif- 
ference between the dereddene d observed o ptical fluxes and 
the appropriate Kurucz model dKurucz 1979) gives the total 
colour excess E(B - V)tot (Table[Tl). Model parameters for our 
sample stars are given in Table [T] and are based on the anal- 
ysis of high-resolution spectra as given by the literature (see 
De Ruvter et al. I d2006h for references). A considerable fraction 
of our sample stars show a photometric variability due to pul- 
sation, so we use only photometric maxima for the SED con- 
struction. 

The errors on the value for E{B - V)toi are calculated using 
a Monte-Carlo simulation on the photometric data. We use an 
error of 0.05 for the photometric measurements in a Gaussian 
distribution. Since we do not know the distances to the sources, 
we adopt a likely luminosity for evolved low-gravity objects of 

= 5000 + 2000 Lq. The distance to IRAS 20056 is likely 
overestimated since visible light from this source probably 
reaches us only by scattering through a nearly edge-on optica lly 
thick disc dMenzies & Whitelocklll98i iGledhill etaDl200 lb . 

We also corrected the infrared spe ctra for reddening by ex- 



tending the average extinction law of ISavage & Mathis ( 
with the theoretical extinction law of ISteenman & Thg d 19891 
199lb and by using the derived total extinction values. This is 
done under the assumption that the extinction is fully due to 
interstellar extinction. Since the total extinction probably con- 
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Table 1. The name, equatorial coordinates or and ^(J2000), spectral type, eff ective temperature T^ff , surface gravity logg and 
metallicity [Fe/H] of ou r sample stars. For the model parameters we refer to iDe Ruvter et al.l j2006l) . Also given is the orbital 
period (see references in lDe Ruvter et al ] l2006l: iGielen et al ] l2007h . The total reddening E(B - V),oi, the energy ratio L^dL, and 
the calculated distance, assuming a luminosity of L» = 5000 ± 2000 Lq. Distances marked with * are upper limits due to the 
aspect angle of the disc. 



N° 


Name 


a (J2000) 


<5 (J2000) 


Tf/-/ 




[Fe/H] 


Pf/'/x'/ 


E(B - V),„, 


LirIL, 


d 






(h m s) 


(° ' ") 


(K) 


(cgs) 




(days) 




(%) 


(kpc) 


1 


EPLyr 


19 18 17.5 


+27 50 38 


7000 


2.0 


-1.5 




0.52±0.01 


3+0 


4.1+0.8 


2 


HD 131356 


14 57 00.7 


-68 50 23 


6000 


1.0 


-0.5 


1490 


0.20±0.01 


50+2 


3.0+0.6 


3 


HD 213985 


22 35 27.5 


-17 15 27 


8250 


1.5 


-1.0 


259 


0.27±0.01 


24+1 


3.1+0.6 


4 


HD 52961 


07 03 39.6 


+ 10 46 13 


6000 


0.5 


-4.8 


1310 


0.06±0.01 


12+1 


2.1+0.4 


5 


IRAS 05208-2035 


05 22 59.4 


-20 32 53 


4000 


0.5 


0.0 


236 


0.00±0.00 


38+2 


3.9+0.8 


6 


IRAS 09060-2807 


09 08 10.1 


-28 19 10 


6500 


1.5 


-0.5 


371 


0.57±0.02 


63+3 


5.4+1.1 


7 


IRAS 09144-4933 


09 16 09.1 


-49 46 06 


5750 


0.5 


-0.5 


1770 


1.99±0.05 


53+5 


2.7+0.6 


8 


IRAS 10174-5704 


10 19 18.1 


-57 19 36 








323 








9 


IRAS 16230-3410 


16 26 20.3 


-34 17 12 


6250 


1.0 


-0.5 




0.56±0.02 


60+3 


6.1 + 1.2 


10 


IRAS 17038-4815 


17 07 36.3 


-48 19 08 


4750 


0.5 


-1.5 


1381 


0.22±0.02 


69+5 


4.5+1.0 


11 


IRAS 17243-4348 


17 27 56.1 


-43 50 48 


6250 


0.5 


0.0 


484 


0.59±0.02 


68+4 


3.8+0.8 


12 


IRAS 19125+0343 


19 15 00.8 


+03 48 41 


7750 


1.0 


-0.5 


517 


1.08±0.02 


52+3 


1.8+0.4 


13 


IRAS 19157-0247 


19 18 22.5 


-02 42 09 


7750 


1.0 


0.0 


120.5 


0.68+0.01 


63+2 


4.2+0.9 


14 


IRAS 20056+1834 


20 07 54.8 


+ 18 42 57 


5850 


0.7 


-0.4 




0.51+0.02 


905+42 


10.9+2.3 


15 


RUCen 


12 09 23.7 


-45 25 35 


6000 


1.5 


-2.0 


1489 


0.55+0.01 


13+1 


2.3+0.5 


16 


SAO 173329 


07 16 08.3 


-23 27 02 


7000 


1.5 


-0.8 


115.9 


0.39+0.01 


36+1 


6.5+1.3 


17 


ST Pup 


06 48 56.4 


-37 16 33 


5750 


0.5 


-1.5 


410 


0.00+0.00 


55+1 


5.7+1.2 


18 


SUGem 


06 14 00.8 


+27 42 12 


5750 


1.125 


-0.7 




0.58+0.02 


111+7 


4.8+1.0* 


19 


SXCen 


12 21 12.6 


-49 12 41 


6000 


1.0 


-1.0 


600 


0.32+0.02 


34+2 


3.8+0.7 


20 


TWCam 


04 20 48.1 


+57 26 26 


4800 


0.0 


-0.5 




0.40+0.02 


42+3 


3.2+0.6 


21 


UY CMa 


06 18 16.4 


-17 02 35 


5500 


1.0 


0.0 




0.00+0.00 


89+3 


9.6+2.0* 



25 - 




20 

Wavelength [/^rn] 

Fig. 1. Comparison between the dereddened (black solid line) 
and reddened (red dotted line) of IRAS 19125. Dereddening 
mainly influences the spectral signature in the 10/im region. 



sists of both an interstellar and a circumstellar component, the 
applied dereddening is thus a maximal correction. 

Dereddening the infrared spectra mainly influences the 
shape of the silicate feature around 9.8/im (see Fig. [T)- Since 
our total extinction values are on average quite small, the effect 
of deredding the infrared spectra is negligible for most stars. 
For stars with a higher value for E{B - V),ot, like IRAS 09144 
and IRAS 19125, one has to be careful interpreting the shape 
of the silicate feature, since the partitioning between interstel- 
lar and circumstellar reddening is not known at this point. 



We compute the energy ratio L/r/L» to determine the 
amount of energy reprocessed by the CE (Table [Ij. 80% of 
our stars have a ratio of 30% and higher, showing that 

the absorption and re-radiation of steUar light by the CE is ex- 
tremely efficient. 

5. General overview 



Figures IbTS] till IrToI show the wide range in observed spectra. 
All spectra are characterised by many distinct spectral struc- 
tures in emission over a smooth continuum. There is quite 
some variety over the sample, with different continuum slopes, 
but for nearly all stars the spectrum is dominated by emission 
bands that can be attributed to silicate dust species. 

In a few stars only, gas phase emission is detected in the 
form of bandhead emission of '^COa and '^COa. These stars 
are HD 52961, EPLyr and IRAS 10174. 

Unique to our sample is the spectrum of EPLyr in the 7- 
20 region. This spectral range is dominated by a strong and 
narrow 11.3 emission peak and a complex plateau with nar- 
row emission features in the 14-18/im region. The strongly 
asymmetric band between 1 jjm and lOyum is real as weU. 
These features are observed in C-rich evolved stars and in the 
carbon-rich component of the ISM and are associated with 
PAH emission and th e associated CH-out-of-plane and C-C-C 
bending modes (e. g. ^ n Kerckhoven et al. I l200(lt iHonv et al 



bOOlHPeeters et aljr2002h . Longward of ISyC/m, the spectrum is 



dominated by silicate emission. EP Lyr, is the only object in our 
sample displaying a mixed chemistry. It is remarkable that in 
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our considerable sample of post-AGB stars, only one possible 
post-carbon star is found! 

5.1. Silicate dust features: Identification 

Amorphous and crystalline silicates are the most commonly 
found dust species in the in terst ellar and circumstellar environ 



ment jMolster et al.lEo02allblB iKemper et alJ 12004 iMin et al 



20071). To describe the glassy and crystalline silicates with an 
olivine and pyroxene stoichiometry we will use the commonly 
used term "amorphous and crystalline olivine and pyroxene". 
Amorphous olivine (Mg2.vFe2(i-A)Si04, where 0<x<l denotes 
the magnesium content) has very prominent broad features 
around 9.8 //m and 18yum (see Fig.|2]), which are easily detected 
in our spectra. These features (also called the lOjum and 20 //m 
features) arise respectively from the Si-O stretching mode and 
the O-Si-O bending mode. For large grains the 9.8//m feature 
gets broader and shifts to redder wavelengths. Amorphous py- 
roxene (MgiFei-cSiOs) shows a 10//m feature similar to that 
of amorphous olivine, but shifted towards shorter wavelengths. 
Also the shape of the 20 fim feature is slightly different. 

The Mg-rich end members of crystalline olivine and pyrox- 
ene, forsterite (Mg2Si04) and enstatite (MgSi03), show strong 
but narrow features at distinct wavelengths around 11.3 - 16.2 
- 19.7 - 23.7 - 28 - 33.6 /zm (see Fig.|2]i, making them easily 
identifiable in our spectra. Forsterite condenses directly from 
the gas-phase at high temperatures (» 1500 K) or it may form 
by thermal annealing of amorphous silicates, diffusing the iron 
out of the lattice-structure. Enstatite can form in the gas-phase 
from a reaction between forsterite and silica (SiOa), or it may 
also form by a similar thermal annealing process as forsterite 



(IBradley et alJ 119831 : iTielens et al.l 119970 . Laboratory experi- 



ments have indicated that silica can be formed when a mor- 
phous silicates anneal to forsterite (e.g. lFabian et al ]|2000). No 
evidence for the presence of silica, with strong features around 
9 and 21 //m, can be found in our spectra. There is also no evi- 
dence for the presence of Fe-rich crystalline silicates, like fay- 
aUte (Fe2Si04). 

Nearly all stars show strong crystalline dust features, both 
at short and long wavelengths, showing that both hot and cool 
crystallines seem to be abundant. 

In our spectra, the amorphous silicate dust seems to peak at 
20 jum, rather than 1 8 fj.m, as would be expected from synthetic 
spectra of amorphous silicates. This could point to the pres- 
ence of Mg-rich amorphous dust, which also shows this shift to 
redder wavelengths (see Fig.|2l). 

5.2. Mean spectra and complexes 

Different complexes at 10, 14, 16, 19, 23 and 33yum can be 
identified in our spectra. In order to study the systematics be- 
tween these complexes in our spectra we compare them to a 
mean spectrum in that region. The mean complex spectra were 
obtained by adding continuum subtracted spectra of sample 
sources with clear spectral structures, using a weighing fac- 
tor proportional to the S/N in that spectral region. The con- 
tinuum was determined by linearly interpolating between the 
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Fig. 2. The mass absorption coefficients [cm^/g] of forsterite 
and enstatite in GRF approximation. Grain sizes of 0. 1 /im (full 
line), 2.0yum (dashed line) and 4.0yum (dotted line) are plotted. 
The third box shows olivine in 0. 1 fim GRF approximation. The 
Mg-rich member (Mg2Si04) is plotted in full line, the standard 
Mg-Fe member (MgFeSi04) is plotted in dot-dashed line. As 
an example we also plotted the SPITZER spectrum of one of 
our sample stars, ST Pup. 



beginning and end of the studied region. The IQfim complex 
is discussed below, other complexes are discussed in the ap- 
pendix (Sect. IB. lb . The mean complex spectra are shown in 
Figs. lB.3l to IB. 71 We also plot the mass absorption coefficients 
of forsterite and enstatite in every complex. 

5.2.1 . The complex (8-13yum) 

In Fig. [3] we plot the continuum divided spectra in the 10 //m 
region, ordered by peak value. We plot 1 + Fy^cJ < Fy^c >, 
where Fy^cs is the continuum subtracted spectrum {Fy - Fy^c) 
and < Fy^c > is the mean of the continuum. The continuum was 
determined by linearly interpolating between 7.5 and 13;um. 
This method preserves the shape of the emission band and al- 
lows a good comparison between the profiles themselves. The 
figure clearly illustrates the very wide variety of the spectral ap- 
pearance of the warm silicates, and because of this large variety 
a mean spectrum was not constructed. We do not find a relation 
between the strength and the shape of the silicate feature. 

In Fig. |4]we show the continuum subtracted flux at 11.3 
and 9.8 /^m versus the peak to continuum ratio of the 10 //m 
silicate feature. The 11.3/9.8/im ratio is a measure for the 
amount of processing that the dust has undergone, where the 
peak/continuum ratio is a measure for the typical grain size. 
We find a high degree of crystallisation (11. 3/9. 8 yum ratio) for 
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Fig. 3. Continuum divided spectra in the lOyum region, ordered 
by peak value. The peak value increases from left to right and 
from bottom to top. The numbering corresponds with the star 
numbers as given in Table[T] For SU Gem and TW Cam we lack 
flux points in the 7.5 jum region to perform a good continuum 
determination. 



all sources, but rather small peak to continuum values in gen- 
eral. EPLyr is a clear outlier with a very strong 11.3/9.8/im 
ratio, due to a very sharp peak at 11.3jum which is due to 
PAH emission. No clear correlation, and thus no evolutionary 
trend, is observed between the shape and the strength of the 
lOyum feature in our sample stars, however. This is in contrast 
with the clear cor relation that is seen in young stellar objects 
( van Boekel et al.H 2003, 2005i), where stars with a strong fea- 
ture have low 1 1 .3/9.8 fim ratios and stars with a weaker silicate 
feature show higher 11.3/9.8yum ratios. Sources that display a 
strong emission feature show rather unprocessed silicate band, 
similar to the ISM profile, with little evidence for crystalline 
grains. Sources with a weak lOyum profile show a broader and 
flatter sihcate feature, dominated by larger grains. 



6. Analysis: Prof ile fitting 

To fit the observed emission features with synthetic spectra, 
mass absorption coefficients for different dust species need 
to be calculated. The conversion from laboratory measured 
optical constants of dust to mass absorption coefficients is 
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Fig. 4. Continuum subtracted flux at 11.3 and 9.8 fim versus the 
peak to continuum ratio of the lOjum silicate feature. SUGem 
and TW Cam are not plotted since for these stars we lack a good 
continuum determination. 



not straightforward and is largely dependent on the adopted 
size, shape, structure an d chemical composition of the dust 
dMin et al.ll20"ol l2005ab . These different dust approximations 
result in very different predicted emission feature profiles, with 
a clear division between the synthetic spectra of homoge- 
neous and irregular grains. If one would assume homogeneous 
spher i cal particles, one could u se Mie theory (lAden & Kerker 
195 1 : Irbon & Ackerman 1981 ). However, cosmic dust grains 
are not perfect spheres, so we have to reside to other meth- 
ods. Examples of irregular dust app roximations are CPE (con- 
tinuous distribution of ellipsoids, BohrenetalJ 1983|). GRF 
(Gaussian random field particles , e.g. iGrynko & Shkuratov 
20031: IShkuratov & Grvnkd l2005h and D HS (distribution of 
hollow spheres, Min et al] 20031 2005 ah grain shapes. The 
DHS shaped particles are further characterised by the frac- 
tion of the total volume occupied by the central vacuum in- 
clusion, /, over the range < / < fmax- The val ue of fmar 



reflec t s the d egree of iiTegularity of the particles (iMin et al 
2003 , 1200 5a'). Cross sections in CDE are computed under 
the assumption that the grains are in the Rayleigh limit (that 
the grains are much smaller than the wavelength of radia- 
tion, thus smaller than 0. 1 i-im), and do not allow the study 
of grain growth. The spectra for the GRF partic les were com - 
puted using the Discrete Dipole Approximation ( Drainell988 ). 
More details on the particle shapes and the c omputations of 
the spectra can be found in iMin et al.l (120071) and Min 2008 
(in preparat ion). Refractory indices f or the materials used ar e 



taken from IServ oin & Pirou (1973); Dorschner et al 
Henning & Stognienko (1996 ); Jaeger et al. Cl998|)~ 



119951); 



Not only grain shape, but also grain size, has a profound 
influence on the mass absorption coefficients. Since a circum- 
stellar disc is the perfect environment for grain growth to occur, 
we also study the effect of different grain sizes on the observed 
emission profiles. Different grain sizes produce emission fea- 
tures at different central wavelengths and with different pro- 
files, as shown in Figs |2] With increasing grain size, emission 
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Fig. 5. Normalised and continuum subtracted emission features 
of the mean of our sample stars, together with mass absorption 
coefficients for different forsterite shape distributions. CDE 
grains are plotted in red, GRF in blue and DHS in green. 0. 1 /jm 
grains are plotted in full lines, 2 /im grains in dashed. 



features will become weaker and eventually disappear, leaving 
mainly a contribution to the thermal continuum. This effect is 
alre ady considerable for grain radii a > 2/im. Previous stud- 
ies ( Bouwman et al. 2001; Honda et"al]|2004i) in the 10/zm re- 
gion have shown that the variety of emission features can be 
described using two typical grain sizes: generally 0. 1 jum to de- 
scribe the grains with radii < 1 fj.m, and l.5-2fj.m for larger 
grains. 

Comparing the observed crystalline emission features with 
calculated synthetic spectra (Fig. |5]l, we find that irregular 
grains are needed to explain the dust profiles. When we look 
at the prominent forsterite features at 16-19-23. 7-33. 6jum, we 
find a good fit is achieved using grains in CDE, GRF or DHS 
(with / = 1.0). None of the tested dust shapes proves an accu- 
rate fit to the, in some stars very strong, forsterite 16 fim feature, 
which seems to be shifted to bluer wavelengths i n the observed 



spectra . This was also seen in the disc spectra of iMolster et al 



(l2002ah . The best fit to this feature is achieved using dust in 
GRF approximation. The fit of the 23.7jum feature is strongly 
improved when using large grains in DHS. The GRF approx- 
imation does not result in a good fit to the 33.6;um feature, 
which is the purest forsterite feature, since at other wavelengths 
the features may be blended with e.g. enstatite emission fea- 
tures. Both CDE and DHS reproduce the shape of this feature, 
but since the CDE approximation does not allow us to study the 
effect of large grains we use GRF and DHS approximations in 
our fitting routine. 

A similar approach to determine the best enstatite or 
amorphous silicate dust approximations is not straightforward. 
Enstatite emission features are mostly blended with forsterite 
features, but in the 14/im region there are some small un- 
blended features. These are ideal features to trace the enstatite 
fraction of the grains. The central wavelengths of these peaks 
are at 13.8, 14.4 and 15.4//m, with the ISAfim feature being 
the most prominent. The sample sources that show clear emis- 



sion in this region all show the 13.8//m feature and in lesser 
degree the 15.4//m feature. The 14.4 jum feature is absent but 
a rather strong unidentified feature appears around 14.8 /im. In 
order to keep the fitting homogeneous we will also use GRF 
and DHS dust approximations to describe the amorphous and 
crystalline pyroxene dust species. 

7. Full spectral fitting 

7.1. Method 

For an exact spectral modelling, full 2D radiative transfer in 
a realistic disc model should be studied. Such models are not 
yet available so as a first approximation we assume the emis- 
sion features to be formed in the thin surface layer of the disc. 
Assuming the flux originates from an optically thin region, 
which is motivated by the fact we see the dust features in emis- 
sion, we can make linear combinations of the absorption pro- 
files to calculate the model spectrum. 

This approach is a first approximation but allows for a gen- 
eral study of trends in dust shapes, grain sizes and processing in 
these discs. Disc models like the one of Dullemond & Dominik 



(120041) . which we used in our pilot study to fit the SED of 
RU Cen and AC Her, do not allow for a detailed spectral mod- 
elling. The code computes the temperature structure and den- 
sity of the disc. The vertical scale height of the disc is com- 
puted by an iteration process, demanding vertical hydrostatic 
equilibrium. Other important processes, like dust settling and 
turbulent mixing, probably also occur in these discs but are not 
included. The code also does not include an independent dust 
species and grain size distribution throughout the disc, making 
a detailed mineralogy study impossible. 

Our model emission profiles are then given by 

F,i~i^aiKdx(Y^PjBATj)) 



where k, is the mass absorption coefficient of dust component i 
and a, gives the fraction of that dust component, B^(Tj) denotes 
the Planck function at temperature Tj and jSj the fraction of dust 
in that given temperature. The temperature of dust grains will 
depend on grain size and grain shape, as well as the distance to 
the central star and the integrated line-of-sight opacity from the 
stellar surface to the grains. This can only be calculated using 
full radiative transfer, so we assume all grains to have the same 
dust temperatures, irrespective of size and shape. 

In this modelling the continuum is given as a sum of Planck 
functions and we take it as another dust component with a con- 
stant mass absorption coefficient. To keep the number of free 
parameters in our fitting routine reasonable, we allow only two 
different dust temperatures and continuum temperatures (be- 
tween 100 and 1000 K), with different ratios. A fit with three 
temperatures was also tested and proved only a minor improve- 
ment. The total number of fit parameters in the model is thus 
15, with contributions of four silicate species, where we use 
two grain sizes, the continuum mass absorption coefficient, two 
dust and continuum temperatures which each their relative frac- 
tions. 
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To fit the spectra we minimise the reduced of the full 
SPITZER spectrum, extended with the IRAS 60 yum flux point 
to constrain the continuum at larger wavelengths, given by 



X 



1 ^ 
~ N - M 2 



■7modei 



with the number of wavelength points Ai, M the number of 
fit parameters and cr,- the absolute eiTor on the observed flux at 
wavelength A,. The errors cr,- represent the statistical noise on 
the spectra. They are generated so that the errors are propor- 
tional to the square-root of the flux and they are scaled in such 
a way that the best fit has a reduced of approximately 1. 

The errors on the fit parameters are calculated using a 
Monte-Carlo simulation. We randomly add Gaussian noise 
with a distribution of width cr,- at each wavelength point. This 
generates 100 synthetic spectra, all consistent with our data, 
on which we perform the same fitting procedure. This results 
in slightly different fit parameters for which we calculate the 
mean, our best-fit value, and the standard deviation. The un- 
certainty on the mass absorption coefficients are not taken into 
account in the;i'^-minimalisation. 

We find that the mass absorption coefficient of small and 
large amorphous grains are quite similar, this degeneracy could 
introduce a large eiTor on the derived fractions for small-large 
amorphous grains. 

7.2. Results 

We tested both the DHS and the GRF dust approximations in 
the fitting routine, where the GRF approximation proved a far 
stronger match. In order to test for the presence of Fe-poor 
amorphous dust, as postulated in Sect. 15.11 we perform the fit- 
ting both with pure Mg-rich amorphous silicates (x = 1, see 
Sect. I5.lt and with the more standard Mg-Fe amorphous sili- 
cate dust {x - 0.5). 

When using 0.1 /urn and 2.0fj.m grain sizes, we find that 
nearly 80% of the dust resides in the 2.0 yum grains. The amount 
of grain growth in crystalline material versus amorphous ma- 
terial is plotted in Fig.|2l Almost all crystalline grains seem to 
be large, with a fraction of large grains in crystalline compo- 
nent higher than 0.7. For the amorphous grains most sources 
have fractions between 0.3 and 0.9. No correlation between 
grain growth in crystalline and amorphous grains is seen. An 
eflicient removal of the smallest grains must have occurred in 
these discs or grain growth was the dominant factor to reduce 
the number of small grains at grain formation. 

Because of this lack of small grains, as a next step we per- 
formed a fitting using grain sizes of 2.0;um and 4.0jum. This 
resulted in a better fit for 17/21 stars, giving slightly better val- 
ues foT x^- For EPLyr and HD 52961 the quality of the fit de- 
creased considerably. This was to be expected, since the strong 
narrow crystalline features in these stars already indicated the 
presence of small grains. 

In general the fit with the jc = 1 .0 Mg-rich amorphous sil- 
icates proved the best, be it only a minor improvement. The 
observed trends in the fit parameters and stellar characteristics 
do not change significantly between the fits when using ;ic = 1 .0 
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Fig. 7. The mass fraction of large grains in the amorphous com- 
ponent versus the mass fraction of large grains in the crystalline 
component, using the fitting with grain sizes of 0. 1 yum and 
2.0 jum. Crystalline grains are almost completely made up from 
large 2.0jum grains. 



or X = 0.5 amorphous dust. Of our sample stars, 6 stars show a 
better fit using the x = 0.5 amorphous silicate dust. 

The best model spectra, overplotted with the observed spec- 
tra are shown in Fig.|6]and Figs. IB.8llB.10l The resulting values 
of the best fit parameters are given in Table lA.fl and Table I A. 2l 
The quality of our fitting is generally very good, 70% of stars 
have ax^ < 5. 

Some trends can be observed (e.g. Fig. |6] and Figs. IB.8I - 
IB.lOl l. As already explained in Sect.|6] the shape of the 33.6 fim 
feature is not well reproduced, although the strength is well 
modelled. Also the 29 fj.m feature seems to be more prominent 
in the observed spectra, than in the model fits. The 19 fim fea- 
ture is slightly overestimated in the model spectra, while the 
neighbouring 23//m feature is slightly underestimated. This 
discrepancy could be due to a data reduction effect, since in 
this region there can be a bad overlap between the short and 
long SPITZER-IRS high-resolution bands. 

The clear outliers in our modelling are EPLyr and 
HD 52961 with very high values of;^^^. These stars have un- 
usual features in their observed spectra, including strong nar- 
row features and CO2 emission lines. A detailed study of these 
outliers will be given in a follow-up article. 

Nearly all sources show the presence of both hot and cool 
dust, both in dust temperature as in continuum temperature. As 
seen in Table lA.ll dust temperatures can differ strongly from 
continuum temperatures. 

7.3. Correlations 

To gain insight in the dust formation process and evolution in 
the circumstellar environment of our sample stars we look for 
trends in the derived fit parameters and correlations with dust 
and stellar parameters. 
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Fig. 6. Best model fits for some of our sample stars, showing the contribution of the different dust species. Best fits for all sample 
stars are given in Figs. IB.8llB.10l The observed spectrum (black curve) is plotted together with the best model fit (red curve) and 
the continuum (black solid line). Forsterite is plotted in dash-dot lines (green) and enstatite in dash-dot-dotted lines (blue). Small 
amorphous grains (2.0/zm) are plotted as dotted lines (magenta) and large amorphous grains (4.0fj.m) as dashed lines (magenta). 



7.3.1 . Mineralogy correlations 

In Fig. [8] the mass fraction in large grains is plotted against 
the mass fraction in crystalline grains. Mass fractions are cal- 
culated as fractions of the total dust mass, excluding the dust 
responsible for the continuum. Note that small grains now in- 
dicate the 2.0 fim grain sizes and large grains the 4.0 fim grains. 

85% of our sources show a mass fraction in large grains 
above 0.5. None of the sources have a mass fraction in large 
grains below 0.25. Most sources have a mass fraction in crys- 
talline grains between 0.1 and 0.6, strongly centred around 0.3. 
There is one clear outlier: IRAS 10174 with a degree of crys- 
tallinity of 0.05. This source shows the more standard ISM pro- 
file, dominated by amorphous grains. No clear correlation can 
be found. 

The fraction of large grains in the crystalline component 
has values between 0.25 and 0.7, while the fraction of large 
grains in the amorphous component has values ranging from 
0.3 till 0.95. Grain growth appears to be more efficient in 
the amorphous dust component. No correlation between grain 
growth in crystalline and amorphous grains is seen however. 

The fraction of enstatite grains in the crystalline component 
lies between 0.1 and 0.5, showing forsterite to be the dominant 



crystalline species. No correlation between the enstatite frac- 
tion and crystallinity can be found. 

In Fig. |9] we plot the continuum to dust luminosity ratio 
of our observed spectra against the mass fraction in large and 
crystalline grains. A correlation can be found between the con- 
tinuum/dust ratio and the mass fraction in large grains. Sources 
with a high fraction of large grains show a high value for 
the continuum/dust ratio. This could be expected since larger 
grains show less prominent emission features and have a larger 
continuum contribution. A high value for the continuum/dust 
ratio also indicates the presence of even larger grains in the 
disc. This indicates that the abundance of these larger grains 
correlate with the abundance of the 4 jum grains, and thus that 
the size distribution continues beyond 4 jum grains. 

7.3.2. Central star correlations 

By comparing the derived fit parameters with stellar charac- 
teristics, like Teff, Lj]f/Lt and the orbital period, we can in- 
vestigate possible evolutionary trends in our sample. However, 
no such correlation between fit parameters and stellar parame- 
ters are found. This was also seen in the SEDs (Sect. |4]i, where 
nearly all stars show a similar SED, irrespective of the central 
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Fig. 8. The mass fraction in large grains plotted against the 
mass fraction in crystalline grains, as derived from our best fit 
parameters. A high degree of crystallinity is found. 




Fig. 9. The continuum to dust ratio of the observed spectra plot- 
ted against the mass fraction on large grains. 

star. Our results do not show an obvious evolutionary trend in 
the mineralogy of our sample stars. 

8. Discussion 

The mineralogy of our sample stars show that the dust is purely 
O-rich. Amorphous and crystalline silicate dust species prevail 
and no features of a C-rich component are found, except the 
PAH emission feature seen in EP Lyr. 

This is remarkable since some of our sample stars are 
thought to have initial masses which would make them evolve 
to carbon stars on the AGB on single-star evolutionary tracks. 
The lack of third dredge-up is also seen in the two objects 
with cl ear CO2 gas emission li nes. EPLyr has a '^C/'^C ratio 
of =s 9 jGonzalez et al. illustratmg that '^C is not en- 

riched during the preceding AGB evolution. The internal chem- 
ical evolution of our sources seems to have been cut short by 
binary interaction processes. 

EP Lyr is the only star in our sample which also has a mixed 
chemistry, and strongly resembles HD 44179. Like HD 44179, 



EPLyr also is strongly depleted dGonzalez et al.lll997bl) . PAH 
emission features dominate the spectrum till 20jum, at larger 
wavelengths crystalline silicates start to dominate. The ob- 
served broad and asymmetric emission feature at 8.2 /vm is very 



simila r to the "class C" objects as discussed in IPeeters et al 



(l20Q2h . The famous "Egg Nebula" also shows this "class C' 
8.2 //m emission feature. EPLyr, together with HD 52961, will 
be studied in detail in a future contribution. 

Some binary post-AGB stars with oxygen-rich discs are 
known that did evolve to carbon stars, like HD 44179, which 
has a large carbon-rich resolved nebula. These stars show in- 
frared spectra which are indicative for mixed chemistries, with 
features of both oxygen-rich and carbon-rich species. The most 
likely scenerio for this is that the formation of the O-rich disc 
antedated the C-rich transition of the central star. Whether other 
of our sample stars will also undergo this evolution is still un- 
known. 

Our full spectral fitting indicates a high degree of dust 
grain processing. The dust seems to consist of considerably 
large grains, with grain sizes larger than 2 /urn. An efficient re- 
moval of small grains must have occuiTed in the discs. The dust 
shape is highly irregular, showing that Mie theory is not ap- 
plicable for the dust in these discs. The spectra of nearly all 
stars show a high degree of crystallinity, where Mg-rich end 
members of olivine and pyroxene silicates dominate. The dust 
condensation sequence of dust in winds of oxygen-rich AGB 
stars predicts the formation of aluminium- or calcium- rich dust 
grains, like co rundum (A1?0^), sp i nel (M gAl?04) or anorthite 
(CaAl2Si208) (' Tielens et all 1997t ICamili2002l) . There is how- 
ever no evidence for the presence of these dust species. 

Most features are well reproduced and only the 14.7 and 
32.5 //m features remain unidentified (Sect. IB. lb . Diopside 
(CaMgSiOs) could be a possible candidate (which has a feature 
around 14.7 fim and a weak feature around 32. 1 fim). However, 
for diopside only mass absorption coefficients for very small 
grains are available, so we cannot include it in the fitting pro- 
cedure in a homogeneous way. 

The dust is highly magnesium-rich, leavi ng a large fraction 



of iron unaccounted for. Previous studies ( Gielen et alJl2007h 



suggest the iron may be locked in the form of metalUc iron. 
Photo spheric depletion i n iron, which we detect i n our sample 



stars jMaas et al.ll2002t I Van Winckel et alj|1998l) . can be un- 



derstood when the iron is locked up in the circumstellar dust 
JWaters etalJl992h . The lack of iron in the detected silicates is 
therefore surprising. If both the crystalline and amorphous sil- 
icates are devoid of iron, t his could mean th at iron is stored in 
metallic iron or iron-oxide ( Sofia et al]|2006() . Metallic iron has 
no distinct features but still makes a significant contribution in 
opacity, especially at shorter wavelengths, making it very hard 
to detect directly. 

In our fitting method, for both dust temperatures equal sil- 
icate fractions are used. This seems to fit features indicative of 
both hot and cool crystalline dust, meaning we do not have a 
strong radial gradient in crystallinity throughout the disc. The 
inner regions of the discs, with temperatures above the anneal- 
ing temperature, are expected to be fully crystalline. The pres- 
ence of a considerable amount of cool crystalline grains implies 
thus that we have strong turbulent mixing in the disc or that the 
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crystalline grains were already abundant at disc formation. It is 
interesting to note that t he very young ci rcumbinary disc in the 
evolved binary SSLep dVerhoelst et al.ll2007i) . which may be 
formed in a similar process as the discs around our post-AGB 
sources, is dominated by small and larg e amorphous grains . 
The crystalline component is very small dSchiitz et al. 2005 ). 
This would indicate that the crystallisation process happens 
during disc evolution and is not already present at formation. 

The evolution of these discs is still unknown. Our anal- 
ysis shows no clear correlation between dust parameters and 
any fundamental parameter (T^/-/^, orbital period) of the cen- 
tral star. In a recent study of IChesneau et al.l (l2007h . a compact 
dusty disc was discovered in "the Ant", a well studied bipolar 
planetary nebula. Interferometric MIDI observations provided 
evidence for a flat, nearly edge-on disc, primarily composed 
of amorphous silicates. This is in contrast with the high crys- 
tallinity observed in the discs around binary post-AGB stars, 
suggesting that the disc in the Ant is relatively young. Whether 
there is a link between this disc and the post-AGB discs re- 
mains unclear. 

Some degeneracy is present in our spectral fits. We per- 
formed the spectral fitting both with pure Mg-rich amorphous 
silicates (x - I) and with the more standard Mg-Fe amorphous 
silicate dust (x = 0.5). Both models often have similar val- 
ues and the best model strongly varies from star to star. Some 
stars also show equally well-fitting models, when using 0.1- 
2.0 yum or 2.0-4.0/im grain sizes. 

8.1. Comparison with young stellar objects 



The mineralogy of the observed spectra shows a striking re- 
semblance to the infrared spectra of young stellar objects, like 
Herbig Ae/Be sta rs (iLisse et al.l|2007|: Ivan Boekel et al.ll2005l) 
or TTauri stars ("Watson et al.' '2007*), and primitive comets 



such as Hale- Bopp (Lisseetal. 2007; Bouw man et al.l 12003 



Min et al.l2005bi) . There also, amorphous silicates and Mg-rich 
crystalline features prevail. The dusty disc is the relic of the star 
formation process so both silicates together with a carbon-rich 
component in the form of PAH emission are often detected, de - 
pending on the disc geometry ( Acke & van den Anckeij|2004 ). 

We compare our findings with the ones discussed in 
van Boekel et al.l (l2005b . They present spectroscopic observa- 
tions of a large sample of Herbig Ae/Be stars in the lOjum re- 
gion. Similar studies on youn g stella r objects have also been 
done by e.g. iBouwman et al.l (1200 lb . One has to be careful 
comparing studies, since our study includes a far larger wave- 
length range, and an exact comparison is thus not possible. 

The degree of crystallinity found in the Herbig Ae/Be 
stars is clearly smaller than for our sample of post-AGB bi- 
naries. The Herbig stars all show a degree of crystallinity 
bel ow 0.35, whereas our st ars have significantly higher val- 
ue s.|van^oekel_etalj(|2005|) use both small (0.1 yum) and large 



(1.5jum) grains to fit the observed emission features and find 
that most sources have a mass fraction in large grains of more 
than 80%. 

For both young and evolved objects a substantial removal of 
the smallest grains has occurred, but seems to be more efficient 



in the discs around the evolved stars. A similar physical process 
might be responsible for the observed grain size distribution 
in both cases, like aggregation of small grains or the removal 
of small grains by radiation pressure. Since several processes 
occur in the formation of the disc it could also be that silicate 
grains in the post-AGB disc are formed at large grain sizes. 
These large grains could be broken up by grain interactions in 
the disc, producing a small fraction of small grains. 



Boekel et al. I (12005') found a clear correlation between 



the crystalhnity and the dominance of enstatite or forsterite in 
these grains. For sources with a high degree of crystallinity 
most crystalline grains appear to be in the form of enstatite, 
while for sources with a low crystallinity, forsterite seems to be 
the dominant species. For all our sample stars forsterite is the 
dominant species, despite the high crystallinity factors. This 
difference could relate to the initial dust species when the disc 
is formed. For discs around YSO, which are formed from the 
ISM, the abundant dust species is likely amorphous olivine. 
Forsterite is expected to be the dominant species formed by 
thermal annealing, while enstatite is expected to be the domi- 
nant species formed by chemical equilibrium processes in most 
of the inner disc. Our discs seem to prefer the formation of 
forsterite, but the initial dust species is not known. The inner- 
most disc regions are hot enough to crystallise dust, but the 
very high degree of crystallinity seems to point to the presence 
of another crystallisation processes, possibly at disc formation. 

9. Conclusions 

We present high-resolution T1MMI2 and SPITZER infrared 
spectra of 21 binary post-AGB stars surrounded by a stable 
Keplerian disc. We summarise our main conclusions: 

- Almost all discs display only O-rich spectral signatures. 
The noticeable exception is EPLyr, which shows a very 
similar spectrum as the central star of the Red Rectangle. 
PAH emission features dominate the spectrum till 20yum, 
at larger wavelengths crystalline silicates start to dominate. 

- Our mineralogy study indicates the dominance of Mg-rich 
amorphous and crystalline silicate dust in the disc. The high 
crystallinity and the large fraction of large grains, as de- 
duced from our full spectral fitting, show strong dust grain 
processing in the discs. 

- The temperature estimates from our fitting routine show 
that a significant fraction of crystalline grains must be cool. 
This shows that radial mixing is efficient is these discs or 
indicate a different thermal history at grain formation. 

- Trend analysis of our fitting parameters show no clear cor- 
relation with stellar characteristics. For the moment it is not 
clear if and how the observed diversity in observed spectra 
relates to specific structural elements of the disc, the star 
and/or the orbits or whether we witness directly an evolu- 
tionary change between different sources. 

To further improve our understanding of these circumbi- 
nary discs, as a next step, we will combine our photo- 
metric and spectroscopic data with interferometric measure- 
ments. Comparing spatial information from the MIDI and 
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AMBER interferometric instruments with a realistic disc 



model dPuUemond & DominikI |2004|) . constrained by photo- 
metric and spectroscopic data, will allow us, not only to study 
the mineralogy, but also the structure of the discs. So far, inter- 
ferometri c measurements for five of our s ample stars have been 
obtained ( Deroo et al ] l2006HDerooll2007h . 



Probing the dust processing in the discs around evolved ob- 
jects proves to be an excellent complement to study physics in 
planet-forming young discs. 
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Appendix A: Tables 



Table A.l. Best fit parameters deduced from our full spectral fitting. Listed the^'^, dust and continuum temperatures and their 
relative fractions. 
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Table A.2. Best fit parameters deduced from our full spectral fitting. The abundances of small (2.0 jum) and large (4.0 //m) grains 

of the various dust species are given as fractions of the total mass, excluding the dust responsible for the continuum emission. 
The last column gives the continuum flux contribution, listed as a percentage of the total integrated flux over the full wavelength 
range. 
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Fig.B.l. The SEDs of our sample stars. The dereddened fluxes (diamonds) are given together with the scaled photospheric 
Kurucz model (solid black line). The coiTesponding TIMMI2 and SPITZER spectrum are overplotted in red. 
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Fig. B.2. See previous caption. 
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B.1. Mean spectra and complexes 

B.2. The 14ium complex (13-1 5 /um) 

In Fig. IB. 31 the 14fj.m complexes are plotted. The mean spec- 
trum has two strong emission peaks, around 13.7 and 14.7/im. 
The profile at 13.7yum seams to blend of two peaks, of which 
one can be identified as enstatite. The small feature in the mean 
spectrum at 14.3 fj.m is also due to enstatite. The feature around 
14.7//m remains unidentified. It seems to be made up of a 
broader feature from 14.3 ;um till 15 fim, with a strong narrow 
peak at 14.7/im. In the sample stars, there is some variation 
in the ratio between the two strong peaks at 13.7 and 14.7 fim, 
meaning they are probably due to different dust species. 
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Fig. B.3. The 14yum complex, continuum subtracted and nor- 
malised. Overplotted in dotted line (red) the mean spectrum. 
The mass absorption coefficients of forsterite and enstatite (in 
CDE approximation) are plotted in green and blue. 



B.3. The IS/j^m complex (15-17 ^m) 

In Fig. IB. 41 the 16/im complexes are plotted. The mean spec- 
trum shows a clear broad feature around 16 fim. The shape and 
strength is very similar to the forsterite feature at 16.2 jim, with 
a contribution of enstatite at 15.3 jum. On top of this broad band 
more narrow features at 15.9/zm and 16.2//m can be discrim- 
inated. The small peaks at 15.4/im and 16.2/zm are due to 
CO2 gas emission, which can only be clearly seen in EPLyr, 
HD 52961 and IRAS 10174. 

Most sample sources show a similar profile as the mean 
spectrum, although the ratio between the two main contribu- 
tors sometimes differs, due to a different enstatite/forsterite ra- 
tio. Stars with a weak 16/zm feature (like TWCam) show a 



very clear separation between the three features, with widths 
of about 0.3 //m. This shows that the mean spectrum probably 
consist of a very broad forsterite feature topped with more nar- 
row features of another dust (or gas?) species. HD 52961 has 
a strong profile that is shifted bluewards in comparison to the 
mean spectrum, with clear CO2 emission lines. 
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Fig. B.4. The Idjim complex, continuum subtracted and nor- 
malised. Overplotted in dotted line (red) the mean spectrum. 
The mass absorption coefficients of forsterite and enstatite (in 
CDE approximation) are plotted in green and blue. 



B.4. The 19 ^m complex (17-21 ^m) 

In Fig. IB. 51 the \9 jjm complexes are plotted. The very broad 
profile is again a blend of forsterite and enstatite. The some- 
times very sharp feature around 19.7 jum is a data reduction 
artefact. Most stars have the same profile shape, except for 
HD 52961, EPLyr and IRAS 10174. The first two are clearly 
outliers in our sample stars, since they have very distinct spec- 
tra, very different from other sample stars. Looking at the 
Spitzer spectrum of IRAS 10174, this star has almost no crys- 
talline features and the observed complex is actually due to 
poor continuum subtraction and normalisation. 

B.5. The 23iJ.m complex (20-27 ^m) 

In Fig. IB. 61 the 23 /im complexes are plotted. The mean pro- 
file is dominated by forsterite and there is an extremely 
good agreement between the mean spectrum and sample stars. 
Remarkable little variation is detected in the whole sample. 



18 



C. Gielen et al.: SPITZER survey of dust grain processing in stable discs around binary post-AGB stars. 




18 19 20 21 18 19 20 21 18 19 20 21 18 19 20 21 18 19 20 21 
Wavelength [/^rm] 



Fig. B.5. The 19/im complex, continuum subtracted and nor- 
malised. Overplotted in dotted line (red) the mean spectrum. 
The mass absorption coefficients of forsterite and enstatite (in 
CDE approximation) are plotted in green and blue. 

B.6. The 33fxm complex (31-37nm) 

InFig. lRTl the 33 jum complexes are plotted. The observed pro- 
file is mainly due to forsterite at 33.6 /zm. At 32.5//m another 
feature can be observed. The ratio 32.5/33.6jum varies slightly 
from source to source, indicating that a different dust species 
is responsible for the feature at 32.5 jum. In ST Pup, the feature 
is clearly broader than in other sample sources. The bump at 
35.7jum is not a reliable result since it sits at the end of the 
spectrum. 
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Fig. B.6. The 23 fim complex, continuum subtracted and nor- 
malised. Overplotted in dotted line (red) the mean spectrum. 
The mass absorption coefficients of forsterite and enstatite (in 
CDE approximation) are plotted in green and blue. 
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Fig. B.7. The 33 //m complex, continuum subtracted and nor- 
malised. Overplotted in dotted line (red) the mean spectrum. 
The mass absorption coefficients of forsterite and enstatite (in 
CDE approximation) are plotted in green and blue. 
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Fig. B.8. Best model fits for our sample stars, showing the contribution of the different dust species. The observed spectrum 
(black curve) is plotted together with the best model fit (red curve) and the continuum (black solid line). Forsterite is plotted in 
dash-dot lines (green) and enstatite in dash-dot-dotted lines (blue). Small amorphous grains (2.0jum) are plotted as dotted lines 
(magenta) and large amorphous grains (4.0/zm) as dashed lines (magenta). 
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Fig. B.9. See previous caption. 



C. Gielen et al.: SPITZER survey of dust grain processing in stable discs around binary post-AGB stars. 

IRAS19157 IRAS20056 




RUCen 




STPup 





15 20 25 30 
X (/Lim) 

SA0 173329 




20 25 
A (/im) 

SUGem 



15 20 25 30 
X (/Ltm) 



40 




15 20 25 30 
X (/Ltm) 



Fig. B.IO. See previous caption. 
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Fig. B. 11. See previous caption. 



